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' In this report, we present the measurements of open charm production 

at mid-rapidity in p+p, d+Au, and Au+Au collisions at RHIC energies. 
I/-) ■ The techniques of direct reconstruction of open charm via its hadronic decay 

f^) ' and indirect measurements via its semileptonic decay are discussed. The 

. beam energy dependence of total charm cross section, electron px spectra, 

and their comparisons to theoretical calculations, including NLO pQCD, 
are presented. The electron spectra in p+p, d+Au, and Au+Au collisions 
at ^/s NN =200 GeV show significant variation. The open charm absolute 
cross section at midrapidity and its centrality dependence are compared to 
those of inclusive hadrons integrated over px > 1.5 GeV/c. 

O ' 1. Introductions 

3 : 

£^ ■ Hadrons with heavy flavor are unique tools for studying the strong in- 

teraction described by Quantum Chromodynamics (QCD). Due to the large 
• ■ charm quark mass, which requires large energies (^3 GeV) for their creation, 

rN \ charm quark production can be evaluated by perturbative QCD (pQCD) 

even at low momentum with the introduction of additional scales related to 
its mass ^H]- Therefore, the theoretical calculations of the charm hadron 
production cross section integrated over momentum space are expected to 
be less affected by non-perturbative processes and hadronization than those 
of the light hadrons [3]. Charm production has been proposed as a sensitive 
measurement of parton(gluon) distribution function in nucleon and the nu- 
clear shadowing effect by systematical studies of p+p, and p+A collisions [3J. 
The reduced energy loss of heavy quarks("deadcone" effect) at momentum 
range 5<pt^,W GeV/c will help us study the energy loss mechanism within 
the partonic medium |5j. A possible enhancement of charmonium (J/ij)) 
production can be present at RHIC energies [HI E] through charm quark 
coalescence. This effect is opposed to the J /if) suppression in a Quark- 
Gluon Plasma(QGP) in the absence of that process jHj- The measurement 
of both open charm yield and charmonia may allow us to quantify the ef- 
fects and study whether charms are in chemical equilibrium with the system. 
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In addition, the heavy flavor transverse momentum distributions and their 
anisotropic flow can be used to study the nature of early thermalization in 
A+A collisions jlUj . 

Identification of charmed hadrons is difficult due to its short lifetime 
(ct(D°) = 124 |im), low production rate and overwhelming combinatoric 
background. Most measurements of the total charm cross section in hadron- 
hadron collisions were performed at low center-of-mass energies ( < 40 GeV) 
in fixed target experiments [111 IT2*] . The measurements at high energy col- 
liders were either at high pj- 13 , with large uncertainty [141 115j or incon- 
sistency [HII- Theoretical predictions for RHIC energy region differ signifi- 
cantly |171 118j . At RHIC, direct measurements of open charm from charm 
hadronic decays and indirect measurements from charm semileptonic decays 
in many beam conditions are possible. We will illustrate the techniques used 
in open charm measurements at RHIC and discuss the results. 

2. Direct open charm reconstruction 

• Large acceptance, Large dataset 

The data used in D direct reconstruction were taken during the 2003 
RHIC run in d+Au collisions at a/s nn = 200 GeV with the Solenoidal 
Tracker at RHIC (STAR) [151120] . A total of 15.7 million minimum 
bias triggered d+Au collision events were used in D° analysis. The 
primary tracking device of the STAR detector is the Time Projection 
Chamber (TPC) which was used to reconstruct the decay of D° — > 
K~ir + (D° — » K + ir~) with a branching ratio of 3.83%. 

• Event Mixing 

The invariant mass spectrum of D°(D°) was obtained by pairing op- 
positely charged kaon and pion in same event with the parent rapidity 
\y\ < 1. The kaon and pion tracks were identified through the ion- 
ization energy loss (dE/dx) in the TPC. The track reconstruction in 
the TPC gives a single track projection resolution of ~ 1 cm around 
the collision vertex and therefore does not allow to resolve the exact 
decay topology. The D° signal in pp < 3 GeV/c and \y\ < 1 after 
mixed-event background subtraction 21 is shown in the left panel of 

Fig.m 

In current STAR analyses, the total charm cross section is largely deter- 
mined by directly reconstructed D(D) at low pt ^3 HO]. In addition, 
D* has been measured at 1.5<pr < 6 GeV/c and L' ± ,D° at pr — 10 
GeV/c |22| . In this presentation, we will focus on total charm cross section 
and its centrality dependence ^3 [23] • 
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Fig. 1. (a) Invariant mass distributions of kaon-pion pairs from d+Au collisions. 
The solid circles depict signal after the mixed-event background subtraction, the 
open circles after subtraction of the residual background by a linear parameter- 
ization, (b) Ionization energy loss (dE/dx) in the TPC vs. particle momentum 
(p) with a TOF cut of |l//3 — 1| < 0.03. Insert: projection on the dE/dx axis for 
particle momenta 1< p <1.5 GeV/c. Figure from [T§| . 



3. Electron spectra from charm semileptonic decay 

There are three techniques used to identify electrons in STAR and 
PHENIX: 

1. Tracking+RICH+EMCAL 

This method is used by PHENIX collaboration to identify electrons |141 
23 . Tracking from drift chamber and pad chamber provides particle 
track with high momentum resolution. A Ring Image Cerenkov Detec- 
tor (RICH) is used as veto counter, and an electromagnetic calorime- 
ter confirms the track's existence and provides a E/p measurement for 
electron selection. 

2. Tracking-r-dE/dx+EMC 

This method is used by STAR collaboration to identify electron with 
p>1.5 GeV/c HEII23. Tb -e STAR TPC not only provides tracking in 
a solenoidal field, but also has good particle identification capabilities 
through ionalization energy loss at resolution of about 8%. Due to 
relativistic rise of dE/dx of electron at high /?7, its separation from 
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7T, K and p are very good at hadron rejection of > 1CP 2 at high px 
. The addition of E/p=1.0 from EMC further reject tt, K and p, and 
hadrons with high dE/dx (nuclei, ghost tracks). Since EMC can be 
used as a trigger detector, the electron spectra from this method can 
be extended to much higher px with high statistics |24| . 

3. Tracking+dE/dx+TOF 

This new technique is developed by STAR collaboration. A prototype 
time-of-fhght system (TOFr) 25 based on the multi-gap resistive plate 
chamber technology was installed in STAR. It covers — 1 < rj < and 
allows particle identification for px < 3.5 GeV/c. In addition to its 
capability of hadron identification , electrons could be identified at 
low momentum (0.2<px < 3 GeV/c) by the combination of velocity 
(/?) from TOFr and dE/dx from TPC measurements. The right panel 
of Fig. ^ demonstrates the clean separation of electrons from hadrons 
using their energy loss {dE/dx) in the TPC after applying a TOFr 
cut of |l//3 — 1| < 0.03. This cut eliminated the hadrons crossing 
the electron dE/dx band. Electrons were required to originate from 
the primary vertex. Hadron contamination was evaluated to be about 
10 — 15% in this selection. 

Gamma conversions 7 — ► e + e~ and ir° —* -ye + e~ Dalitz decays are the 
dominant photonic sources of electron background. There are again several 
methods used to study the background electron sources: 

1. Cocktail Modeling 

This method was adapted by PHENIX Collaboration in their first 
publication 14 j. The cocktail is a detailed detector simulation of 
electrons/positrons from 7 conversion and Dalitz decays from 7r°,r] 
and other photonic sources. 

2. Converter 

This method was adapted by PHENIX Collaboration in their latest 
report Two datasets with and without a converter with known 

thickness were taken and analyzed. The two measurements allows 
to untangle the sources associated with photon conversion and those 
without. Species-dependent Dalitz fraction per 7 and contribution 
from other photonic sources (j],u), p,(j> and K) are evaluated by de- 
tailed detector simulations. 

3. Invariant Mass Reconstruction 

The e + e~ pairs from photon conversions and Dalitz decays are present 
mainly at small pair invariant mass and/or small opening angle |191 
I20j . Due to the large coverage of the TPC, the efficiency of finding 
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Fig. 2. Non-photonic electron px distributions from p+p, p+p, d+Au,and Au+Au 
at different beam energies from ISR to Tevetron. 

pairs is very high for such processes. To measure the background, 
the invariant mass and opening angle of the pairs were constructed 
by first selecting an electron (positron) in TOFr and then matching it 
with every positron (electron) candidate reconstructed in the TPC 
without additional requirement of a secondary vertex at the conver- 
sion point. In this method, both inclusive electron spectra and back- 
ground sources are taken from same dataset, and the reconstructed 
background spectrum is not sensitive to the detailed knowledge of the 
conversion and Dalitz decay. About 95% of electrons from sources 
other than charm semileptonic decays have been measured with this 
method, while the remaining fraction (< 5%) from decays of rj, u, p, (ft 
and K was determined from detailed detector simulations. 



4. Discussions 

Fig. |2] shows the energy dependence of the electron px spectra from 
ISR to Tevetron. Direct open charm measurements (STAR, CDF) are con- 
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Fig. 3. Total cc cross section per nucleon-nucleon collision vs. the collision energy 
(in A/A- 



verted to electron spectra. Several measurements at RHIC have shown to 
be consistent between (dE/dx+EMC), direct D, (dE/dx+TOF) [2011221121] 
at STAR, and to some extent, the measurements in p+p collisions between 
STAR and PHENIX 27 . Fig. [2J shows significant dispersion of electron 
distribution at ^/s NN =200 GeV at RHIC from Au+Au, d+Au to p+p at 
Pt — 2 GeV/c where statistics are still high. The variation in spectra is 
possibly due to the interactions between the charm quark (hadron) and the 
medium [2Hj. The deadcone effect would be most effective at 5^pr~10 
GeV/c and not yet accessible at the px range (~ 2 GeV/c) of electrons. 
The beam energy dependence of the cross section is shown in Fig. |SJ The 
PHENIX results are derived from electron spectrum in Au+Au minbias 
collisions using the electron spectrum shape from PYTHIA6.205 }2rS| . The 
STAR result is from a combined fit of D° and electron spectra in d+Au 
minbias collisions. At y/s ~ 52 — 63 GeV, the available measurements are 
inconclusive due to the inconsistency between different measurements 
and are omitted. At ^/s^ = 200 GeV, both PYTHIA and NLO pQCD 
calculations underpredict the total charm cross section |17l I29j . This is 
evident from STAR data but less pronounced in PHENIX results. There 
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are indications that a large charm production cross section (o^ — 2 — 3 
mb) at Js NN ~ 300 GeV is essential to explain cosmic ray data [30]. The 
centrality dependence of charm da/dy per binary nucleon-nucleon collision 
at midrapidity shows possible dependence on centrality with lower produc- 
tion cross section in central collisions as in Fig. In the same figure, the 
cross sections per binary nucleon-nucleon collisions of inclusive hadrons |31j 
(da h ±/dr]/3(fj,b)), which were integrated over px > 1.5 GeV/c (close to the 
mass of charm hadron mp) and scaled down by a factor of 3 for three light 
quark flavors, show same magnitude of absolute cross section and a similar 
trend. This may suggest that particle production rate is not sensitive to 
flavor when the momentum transfer is above the production threshold. 
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Fig. 4. Charm differential cross section per binary nucleon-nucleon collision 
(d<j^ e N j dy) at midrapidity as function of number of participant micleons at RHIC 
y / s7^7=200 GeV. Also plotted as open stars are cross section per binary nucleon- 
nucleon collisions {du^± /dn of inclusive hadrons integrated over px > 1.5 
GeV/c (~ mo) and scaled down by a factor of three. 



A direct reconstruction of open charm spectra in Au+Au collisions from 
0^Pt~10 GeV/c in future runs, when compared with those from d+Au 
data |191 122 j . will enable us to study the thermalization at low px and 
deadcone effect at high px ■ The elliptic flow V2 measurements of charm 
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flow with higher statistics at the same px range will further strengthen the 
case, since early thermalization will result in large v 2 at low px , and absence 
of large energy loss due to deadcone effect may result in smaller v-i at high 
PT when compared to those of light hadrons. 

5. Summary 

In summary, charm cross section and transverse momentum distribu- 
tions from p+p, d+Au and Au+Au collisions have been measured by the 
STAR and PHENIX collaborations at RHIC. The NLO calculation under- 
predicts the total cross section in this energy range. A possible centrality 
dependence of charm cross section can explain the marginally different cross 
sections observed by STAR in d+Au collisions and PHENIX in Au+Au col- 
lisions. We look forward to more results using charm as probe to the hot 
and dense medium created at RHIC from run4 and future runs. 
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